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Abstract —The free-space optical (FSO) communications can 
achieve high capacity with huge unlicensed optical spectrum and 
low operational costs. The corresponding performance analysis 
of FSO systems over turbulence channels is very limited, espe¬ 
cially when using multiple apertures at both transmitter and 
receiver sides. This paper aim to provide the ergodic capacity 
characterization of multiple-input multiple-output (MIMO) FSO 
systems over atmospheric turbulence-induced fading channels. 
The fluctuations of the irradiance of optical channels distorted by 
atmospheric conditions is usually described by a gamma-gamma 
(FF) distribution, and the distribution of the sum of FF random 
variables (RVs) is required to model the MIMO optical links. We 
use an a-/i distribution to efficiently approximate the probability 
density function (PDF) of the sum of independent and identical 
distributed FF RVs through moment-based estimators. Further¬ 
more, the PDF of the sum of independent, but not necessarily 
identically distributed FF RVs can be efficiently approximated by 
a finite weighted sum of PDFs of FF distributions. Based on these 
reliable approximations, novel and precise analytical expressions 
for the ergodic capacity of MIMO FSO systems are derived. 
Additionally, we deduce the asymptotic simple expressions in high 
signal-to-noise ratio regimes, which provide useful insights into 
the Impact of the system parameters on the ergodic capacity. 
Finally, our proposed results are validated via Monte-Carlo 
simulations. 

Index Terms —Free-space optical communications, achievable 
rate, atmospheric turbulence, gamma-gamma distribution. 

I. Introduction 

Free-space optical (FSO) communication is one of the 
most promising technologies for indoor and outdoor wireless 
applications due to its merits of free license, effective cost 
and high bandwidth. The research of FSO systems has been 
taken around for the last three decades and is currently 
attracting more attention as the demand for high data rate 
continues to increase Q, m. The typical applications of FSO 
systems include “last mile” access, indoor positioning, disaster 
recovery, military applications, underwater system, device-to- 
device communications, and video transmission, etc. 11-®. 

However, the FSO communication still faces many chal¬ 
lenges. One of the main challenges arises from the atmospheric 
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turbulence, which is induced by the fluctuations in the atmo¬ 
sphere because of inhomogeneities in pressure and temperature 
changes, especially for over link distances of 1 km and above 
0. The performance of FSO systems is highly vulnerable to 
the atmospheric channel turbulence between the transmitter 
and the receiver, and the parameters of the channel such as 
the length, the optical wavelength, and turbulence. Therefore, 
accurate modeling for the distribution of the turbulence- 
induced fading is significant to assess the performance of 
FSO communications. To describe weak or strong atmospheric 
turbulence fading, the gamma-gamma (FF) distribution has 
been proved to provide good agreement between theoretical 
and experimental data Q-iISl. 

The performance of FSO systems over FF fading chan¬ 
nels has been widely investigated in the literature lfT0ll - llT4| . 
These research have proved that effective fading-mitigation 
techniques are required to satisfy the typical bit error rate 
(BER) and capacity targets for FSO applications at the 
range of practical signal-to-noise ratios (SNRs). One of the 
promising approaches to mitigate the degrading effects of 
atmospheric turbulence is the use of multiple-input multiple- 
output (MIMO) technique, which has been widely used in 
wireless radio systems. By deploying multiple apertures at the 
transmitter and/or the receiver, the FSO system performance 
can be significantly enhanced. Therefore, many researchers 
have focused on the performance analysis of MIMO FSO 
systems over turbulence channels ifTSll - ll^ . In ifTSll - lflTll . the 
authors have studied the BER, diversity gain and combining 
gain of MIMO ESO systems with equal gain combining (EGC) 
or maximal ratio combining (MRC) at the receiver. The outage 
probability of ESO communication systems over FF fading 
channels with spatial diversity has been provided in lfT8ll - ll20| . 

The performance of MIMO FSO systems over FF fading 
channels has been extensively investigated in terms of outage 
probability and error rate, however, there have been few studies 
on the theoretical ergodic capacity. This is due to the fact that 
the integral for ergodic capacity is mathematically difficult 
to be solved. Only very recently, the ergodic capacity of 
MIMO FSO communications using multiple partially coherent 
beams propagation through strong turbulence channels has 
been investigated in ET\ by using a single FF approximation. 
However, the approximation results presented in need 
an adjustment parameter to obtain a sufficient approximation 
accuracy. In this paper, we use a more accurate approximation 
method to evaluate the capacity of MIMO FSO systems. 

Capitalizing on the aforementioned observations, in this 
work we provide the cumbersome statistical ergodic capacity 
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analysis of MIMO FSO systems with EGC over turbulence 
channels, which are modeled as independent and identically 
distributed (i.i.d.) and independent, but not necessarily iden¬ 
tically distributed (i.n.i.d.) FF distributions, respectively. In 
particular, the contributions of this paper are summarized as; 

• Accurate analytical expressions for the ergodic capacity 
of MIMO FSO systems with EGC over i.i.d. FF fading 
channels have been derived. These expressions are given 
in terms of either Meijer’s G- or Fox’s iT-functions. 
The Meijer’s G-function can be easily evaluated and 
efficiently programmed in most standard software pack¬ 
ages (e.g., MAPLE, MATHEMATICA), while numer¬ 
ically efficient methods can be used to evaluate the 
Eox’s 7J-function fT2\ . ll^ . The presented technique of 
ergodic capacity analysis is founded on approximating 
the probability density function (PDE) of the SNR at the 
receiver by using the a-fj, distribution ll24ll . Compared 
with the approach employing in 1211 . the a-/j, approxi¬ 
mation of the sum of multiple FF random variables (RVs) 
is more accurate and avoid a correcting factor to obtain 
a sufficient approximation accuracy. 

■ Based on the nested finite weighted PDF of the sum 
of i.n.i.d. FF RVs presented in ||25]| . we also derive a 
novel and analytical expression for the ergodic capacity 
of MIMO FSO systems for the case of i.n.i.d. FF fading 
scenario. 

• To further provide the useful yet simple insights into the 
impact of system and channel parameters on the ergodic 
capacity, we present the asymptotic ergodic capacity 
bounds of MIMO FSO systems over i.i.d. and i.n.i.d. 
FF turbulence channels in the high-SNR regime. For 
example, it shows that the ergodic capacity increases 
as the number and diameter of apertures gets larger. 
Interestingly, it is not a proportional relationship between 
the ergodic capacity and link distance as expected. 

The remainder of the paper is structured as follows: Section 
El describes the system model and the FF distribution of the 
turbulence fading. In Section [Till we present the approximated 
PDFs of the sum of multiple i.i.d. and i.n.i.d. FF RVs, and 
derive exact and high-SNR approximation expressions for 
the ergodic capacity. Numerical and Monte-Carlo simulation 
results are provided in Section IIVI and Section |V] concludes 
this paper with a summary of the main results. 

Notation: We use upper and lower case boldface to denote 
matrices and vectors. The expectation and variance is given 
by E(-) and Var(-), respectively. The real part of a number 
is given by 5P(-). The matrix determinant is given by det(-), 
while the Hermitian operation of a matrix X is defined as 
. The set of positive integer numbers is expressed via Z+. 
Finally, let Xij denote the (z, j)-th element of the matrix X. 

II. FSO System Model 

We consider a typical outdoor MIMO FSO system with M 
transmit lasers and N receive apertures with same hardware 
over atmospheric turbulence channels. Similar to the related 
literature oa, Eii-Ea, we assume that the high-energy FSO 
system adopts the commonly used intensity modulation/direct 


detection(IM/DD) of on-off keying (OOK) signals using in¬ 
terleaving coding. In this case, the MIMO FSO system’s 
performance is limited by background radiation and thermal 
noise, which can be modeled as i.i.d. additive white Gaussian 
noise (AWGN) as a accurate approximation of the Poisson 
photon-counting detection model ifTSl . ifTa . ||2TI . Il27l . Then, 
the received signal at the nth receive aperture can be expressed 
as 

Un — ( 1 ) 

where s G { 0 , 1 } denotes the transmitted signal, Pavg is 
the transmitted average power, p represents the optical-to- 
electrical conversion coefficient, and Vn is the AWGN with 
zero mean and variance of = No/2. Moreover, = 
Sm=i denotes the sum of M fading gains at the nth 

receive aperture, where is the irradiance from the mth 
transmit laser to the nth receive aperture. Since the typical 
FSO system provides high data rate (about 10® bits/s), the time 
scales of these fading processes Imn are far larger than the bit 
interval (about 10“® s). The ergodic fading channel, in which 
each transmitted codeword experiences the entire statistics 
of the fading process, is commonly assumed to analyze the 
theoretical performance of FSO systems Qoi, im. Moreover, 
we assume the fading gains are independent RVs. This can 
be easily achieved by placing the apertures a few centimeters 
apart, due to the fact that the coherence length of the irradiance 
fluctuations is of the order of centimeters insi, m- 


A. Electrical SNR 


We normalized the average gains Imn — ^{Imn) at each 
branch to unit. It has been shown in m that in turbu¬ 
lence fading channels, the performance gains provided by the 
complicated MRC scheme are only small to moderate when 
compared with the simple EGC scheme. Therefore, EGC is 
more attractive in ESO systems due to its remarkable lower 
implementation complexity. When the EGC method is used, 
the received sum output signal is 


y = 


'/Pys 

MN 


I + v, 


( 2 ) 


N M 

where I — Imn, and v is the effective Gaussian noise 

n=l m—1 

at the receiver. We normalized the total transmit power P = 
MP^vg = 1- Note that the division by MTV in (|2]l is used to 
ensure that the total transmit power and the sum of the receiver 
aperture areas are the same as that of a system with single 
transmit/receive antenna, respectively. Then, the instantaneous 
electrical SNR of the sum output signal at the receiver can be 
calculated as 


7 = 


{yif 

{MNfNo 


7o/^ 


(3) 


where 70 = with p = I/Nq denoting the transmit 

SNR, and the average electrical SNR can be correspondingly 
denoted by 7 = 70 /®. 
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B. rr channel model 


by using ||29] Eq. (6.561.16)] as 


Adopting suitable model for the turbulence fading channel 
is of key importance to assess the performance of FSO 
systems. It has been proved in m, I?) that the EE distribution 
can well characterize the fading gains Imn for turbulence 
scenarios from weak to strong. Based on a doubly stochastic 
theory of scintillation, the small-scale irradiance fluctuations 
are modulated by large-scale irradiance fluctuations of the 
propagating wave. Then, the EE distribution can be derived 
from the product of two independent Gamma distributions and 
the PDF of Imn is given by ll25l 


//„„ (/mn) = 


2{aby 


E (a) E (6) Im 
X Ka-b 1 2 



(4) 


where a and b are the shaping parameters of small-scale and 
large-scale eddies of the scattering environment, E (•) denotes 
the Gamma function ||29] Eq. 8.310.1], and K„ (•) is related 
to the modified Bessel function of the second kind with order 
V ll^ Eq. 9.6.1]. In this paper, we assume the spherical wave 
propagation, which is the exact model for line-of-sight MIMO 
systems OTI and widely adopted in the literature ifTOll . ifT^ . 
Ha, im, ES). The effective numbers a and b are related to 
the atmospheric conditions and are respectively given by m 


/ 


exp 


0.49cj| 




7/6 


b = 


exp 


(l-f0.18d2-f0.56cr. 
0.51cr|(l-f 0.69crf/®) 

(l-f 0.9d2-f 0.62d2cr2^/^)^^'^ 


- 1 


- 1 


(5) 


( 6 ) 


where erf = 0 .492C2fc7/6j^ii/6 is the Rytov variance and 
d = ‘\JkD'^/AL with L being the distance between transmitter 
and receiver, k = 27r/A denotes the optical wavenumber with 
A being the operational wavelength, D is the aperture diameter 
of the receiver, and is the altitude-dependent index of 
the refractive structure parameter determining the turbulence 
strength. Furthermore, we assume remains constant for rel¬ 
atively long transmit bits interval and varies from 10“^^m“2/3 
to 10“^^m“2/3 for weak to strong turbulence cases ESI. 

According to the experimental data, we can properly choose 
the values of a and 6 in dUi to provide a good approximation of 
the irradiance fluctuation PDF. Moreover, a and b are linked to 
the scintillation index, which is used to describe the strength 
of atmospheric fading and given by ll^ Eq. (4)] 


S.I. = 


E(J^) . _ 1 I 1 

[E(/)]2 a b 


1 

ab 


(7) 


From (|7]l, it is clear to see that the scintillation index depends 
on the values of a and b. Furthermore, the turbulence becomes 
weaker for decreasing scintillation index as a and/or b increase, 
and gets stronger for increasing scintillation index with smaller 
values of a and/or b. The qth moment of Imn can be obtained 




T {a + q)r {b + q) Iab\ 

r(a)r(6) 


( 8 ) 


III. Ergodic Capacity Analysis 


Considering the number of transmit apertures is less than 
the number of receive apertures, e.g., M < N, the ergodic 
capacity of MIMO FSO system in bits/s/Hz is given by 



where Im denotes an MxM identity matrix, and H is an Nx 
M channel matrix with elements being Inm- In this section, we 
present some statistical features of the sum of L = MN i.i.d. 

L 

and i.n.i.d. EE RVs where S' = ^ /;. The approximation 

methods for the PDF of the sum of EE RVs are introducecQ. 
These results are useful for the ergodic capacity analysis of 
MIMO FSO systems over turbulence channels. 


A. I.i.d. EE Fading Channels 


Recently, the a-fi distribution has been proposed to approx¬ 
imate the sum of multiple i.i.d. EE RVs |[35l . Ii3^ . This is 
because the a-^L distribution is a more general distribution, 
which includes the Gamma distribution as a special case. We 
use R to denote an a-/i RV. Then, the PDF expression of a-/i 
distribution is given by li2^ 




1 

f“i^E (/r) 


exp 



( 10 ) 


where a > 0, /r = E^ (77“) /Var(i?“) is the inverse of the 
normalized variance of i?“, f = [E is a a-root mean. 

With the knowledge of the first, the second and the fourth 
moment of R, the moment-matching method can be used to 
approximate the parameters a and p, ll24l . Therefore, we need 
to derive the parameters a, p and f by solving the following 
nonlinear functions: 


(/) 

E(/2)-E2 (/) 

E^ in 

E(/4)-E2 (/2) 


E(p)E(p + 2/a)-E2 (p+l/a)’ 
E2 (p-f 2/a) 

E (p) E (p + 4/a) - E2 (p + 2/a) ’ 
pV°E(p)E(/) 

E(p + l/a) ’ ^ 


where the moments in (fTTI) can be evaluated by using ([8]) and 
the multinomial identity as lf35] Eq. (10)] 


31=0 31 - 1=0 


q 


Jl-2 

jl-1 


X E 




( 12 ) 


^The PDF analysis can also be studied by the method of generalized power 
series representation o, which induces an infinite series PDF expression. 
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Note that an analytical solution to the equations (fTTT i is very 
difficult to be obtained, so we use numerical methods instead, 
such as the fsolve function of Matlab and Maple. 


tions can be expressed in terms of the Fox’s i?-function as 


In (1 + x) = 


( 1 , 1 ), ( 1 , 1 ) 
( 1 , 1 ), ( 0 , 1 ) 


( 20 ) 


Based on (|9]), the ergodic capacity of MIMO FSO systems 
with EGC can be expressed as 


C = 


pCO 

/ log 2 (1 + 7 ) / ( 7 ) fi 7 . (13) 

Jo 


Using a simple power transformation of Q, we have 


/(7) =/fi(7)-r- (14) 

d'f 

Recall (doll, the PDF of the instantaneous electrical SNR 7 at 
the receiver side can be given by 


= H, 


1,0 


0,1 


( 0 , 1 ) 


( 21 ) 


Substituting (l2lT i and (l20l i into (fThl) . we get the Mellin-Barnes 
integral ll^ of the product of two Fox’s iF-functions. By 
involving Eq. (2.25.1.1)], the ergodic capacity of MIMO 
ESO system with EGC scheme for arbitrary values of a and 
fj, is given by 


Ci.i.d. = 


1 

F (/r) In 2 


H- 


1,3 

3,2 


'1^10 

( 1 , 1 ), ( 1 , 1 ), ( 1 -/ 4 , 2 /a) 

^2/a 

( 1 , 1 ), ( 0 , 1 ) 


( 22 ) 


/i.i.d. ( 7 ) 


ap^7 2 


27o" (P) 


■ exp 



(15) 


where I is the a-root mean value of the approximation. Then, 
di can be rewritten as 

27 “'"/^/“/^r(p)ln 2 

X J 7°^/^-i In (1 + 7 ) exp ^ ^ 7 . (16) 


To evaluate (fThl l. we can express the logarithmic and expo¬ 
nential functions as Meijer’s G-functions by IIJTI Eq. (11)] 


In (1 -h a;) = 


1,1 

1,0 


-X _ (^1;0 
, ^ — '-^0,1 


0 


(17) 


Note that (|22| is very compact in the form of Eox’s H- 
functions and will be verified through Monte-Carlo simula¬ 
tions in Section |IV] We can use the numerically efficient 
methods presented in ll22ll . 12^ to evaluate the Eox’s H- 
function. The path of the integration is running from w —too to 
w-l-too, which depend on all poles of gamma functions in the 
numerator iQi. The ergodic capacity expression of SISO ESO 
links has been given in IfTOl Eq. (16)] by using the PDE of FF 
distribution directly. Flowever, the a-p approximation method 
is utilized in our work to study the MIMO FSO system. It 
should be pointed out that although (| 22 | cannot intuitively 
reduce to the existing expression for the SISO case (e.g., 
mni Eq. (16)]), their mathematical calculating results are quite 
close to each other, which will be verified by the simulation 
results in Eig. 1 and 2 (M = N = 1 for the SISO link). 


Upon using IIJTI Eq. (21)], the analytical expression of ergodic 
capacity is given by 


Ci.i.d. = 


a/r' 




2(7o" /“/^F (fi) {2 tt 


l+kl2-3j2 


In 2 


X G 


k+2l,l 

2l,k+2l 


[( 

A0,-^),A(U-T) 

\ki^-fS ) 

Aik,0),A{l,^),A{l,^) 


(18) 


B. I.n.i.d. FF Fading Channels 

The ergodic capacity for i.i.d. FF fading channels has been 
derived in (i 22 l l. now we move on to the ergodic capacity 
analysis for the i.n.i.d. FF fading channels. It is reasonable 
to assume i.n.i.d. fading when the transmit or receive aper¬ 
tures are placed far from the others. The approximative PDF 
expression for the sum of multiple i.n.i.d. FF RVs has been 
presented in li^ by neglecting the cross terms as follows 


where A (e, r) = • • • , with r is an arbitrary 

real value and e is a positive integer. Moreover, l/k = a/2, 
where I and k are both positive integers. For example, if a = 
0.8, we can set / = 2 and fc = 5. Furthermore, I = a and k = 2 
are for the special case of a G IP. Note that the calculation 
of (fTsT i is efficient for special values of a. However, for large 
values of I and k, it is very tedious to use (iTSi) . Therefore, 
another method is adopted to solve this problem as follows. 


We recall the well-known translation from the Meijer’s G- 
function to the Fox’s iF-function as Eq. (8.3.2.21)] 


^p,q 

X 

[a^, 1] 

\h 11 

_ ^m,n 

— ^p,q 

X 



[Gpf IJ 




(19) 


With the help of (fTTl i. the logarithmic and exponential func- 


5 = = (23) 

1=1 1=1 \l=l / \l=l J 

where xi and yi are Gamma RVs with parameters {ki,l/ki) 
and respectively, and // = xiyi is a decom¬ 

position from a FF RV to two independent Gamma RVs. 
Following the well-known property for the product of two 
random variables 13^ . the PDF of /; can be given in the 
form of (in with the shape parameters ki = a, mi = b, and 
the mean power U; = Imn- In most practical MIMO FSO 
systems, the transmit lasers and receive apertures are placed 
only a few centimeters apart. Then the large-scale irradiance 
fluctuation changes relatively slowly compared with the small- 
scale effect. It is reasonable to assume that the FF RVs have 
one shaping parameter in a column, i.e., ki = k. Based on 
(I 23 I 1 . the PDF for the sum of L i.n.i.d. FF RVs is shown in 
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Il25l Eq. (27)] as 

L rtii 

fs{s) = '^'^WL 

i=i j=i 


r(Lfc)r(j)^ 


Lk+j 


~^Lk—j 


Lkrrij 


( 


where the weights can be easily computed by the recursive 
formula according to ||25] Eq. (24)] 

wl {i, rrii - t, 


^ L t 

= 7 E E 


n^rriq 
rUn [ 1- ——- 
\l„mi 


9=1.94* i=l 
X WL {i,mi -t + j, 

where t = 1, - ■ ■ , — 1 , and 


(24) 


Using IIJTI Eq. (21)], the integral in (|26| can be expressed as 


J ln(l + 7 ) 7 " d 7 


7 U 2 2 


7 


1,1 

1,0 


^ 2,0 

'-^ 0,2 


^\/7 


V_ —V 
2 ’ 2 


_ 1 ^6,1 
47r^2.6 


E 

16 


—U, 1 — It 

v_ v -\-2 y_ v — 2 

4’ 4 ’ 4’ 4 ’ 


d7 

(28) 


where we have used the notations of it = [Lk + j)/4, v = 
Lk — j, z = Lfcmi/(fli^/To), and the identity ll^ Eq. 9.34.3] 


K,(x) 


= 


. 2,0 

0,2 




— 

T 

V 

2 

V 

2 . 


(29) 


By substituting (l28T l into (l26l l. the ergodic capacity of MIMO 
ESO systems over i.n.i.d. EE fading channels is expressed as 
( l30l ) at the bottom of this page. 




L 


n 

1=1 j'4* 



Eor MIMO ESO systems with EGC receiver, we can ap¬ 
proximate the PDE of the instantaneous electrical SNR 7 when 
one of the shaping parameters remains the same with a nested 
finite weighted sum of EE PDEs as follow: 


L rrii 


/i.n.i.d. ( 7 ) = E E^^-^ {w/}f=l, {U/}f=i) 


i=l j=l 

Lk+j Lk+j 

{Lkmi) = ^7 2 

E(j)E(Lfc)(aV7^) ^ 


-Khk-I 


I Lkrrii'^ 

Uj7o 


(25) 


Note that the accuracy of (i25l l depends on the combinations 
of the parameters k, mi, and O;. This constraint has been 
explicitly investigated by employing Kolmogorov-Smirnov 
goodness-of-fit test in Section III. C of ll25l . 


Eor i.n.i.d. EE turbulence channels, the ergodic capacity in 
([I3]l can be written as 

L rrii 

C....a. = EE CijWL 

i=i j=i 


ln(l+7)7^-'KLfe-, 


2 / Lkm,+j 
. \ 


d 7 , (26) 


where Cjj is a constant dependent on z, j, and given by 


A 

C/ j — 


Lkmj 


E(LA:)E(j)ln2 


(27) 


C. High-SNR Analysis 


In order to get more insights into the impact of system and 
channel parameters on the ergodic capacity, we investigate the 
asymptotic ergodic capacity of ESO MIMO systems in the 
high-SNR regime. Eirst, we consider the i.i.d. case. By taking 
7 large in (fThI ). the ergodic capacity at high SNRs becomes 


C 


00 _ 

i.i.d. 


1 /*00 n iLA —1 

1 / , / t 

r+ - 

111 2 Jo 27 o"/“/^E(p) 




(31) 


With the help of ED Eq. (2.6.21.2)] 

exp (—px®) \nxdx 

= e-^p-G^V (C/0) [ 1 /: {CIO) - Inp], (32) 



(ED can be further calculated as 


C OO _ 

i.i.d. — 


;ln2 - 


a 


(p) — In p -I- a In / -f — In 70 


(33) 


where ip (•) is the Euler’s digamma function ll^ Eq. 
(8.360.1)]. It is obvious that the high-SNR ergodic capacity 
increases with the optical-electrical conversion coefficient p. 
Note that similar observations were also made in ED, ED . 


Now we consider the high-SNR ergodic capacity of MIMO 
ESO systems over i.n.i.d. EE fading channels. Using the 
similar method aforementioned, we approximate ln(l-|- 7 ) with 
ln( 7 ) by taking 7 large in the i.n.i.d. case, so the integral in 


Lk+j 


L rrii 

EE 

i=i j=i 


1 

' Lkmi 


47rE(Lfc)E(j) In 2 

\^ilo 

xGlll 

( 1 

/ Lkmi\^ 


Lk+j 1 Lk+j 

4 5-*- 4 

1 I 670 

[ n, ) 

Lk-j 

4 ’ 

Lk—j+2 Lk—j Lk—j—2 Lk-\-j Lk+j 

4 ’ 4 ’ 4 ’ 4 ’ 4 


( 30 ) 
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(|28] t can be expressed 


111 ( 7 ) 7 " 


X r 


Au — V 




Au + V 


<17 = 

/4 
■f 


z- 

4u — V 


4u + V 

2 

-ln(z) 


(34) 


where we have used the identity in Il43l Eq. (2.16.20.1)]. 
Substituting 01 into 01 and after some basic algebraic 
manipulations, the analytical expression for the high-SNR 
capacity can be hnally derived as 


L rrii 


.d. = 

i=i i=i 

X [?/) {Lk) + i’ (j) - in {Lkrrii) + in . (35) 


The above expression is intuitive, since it indicates that at 
high SNRs, the larger fading power Vli, transmit SNR p, 
and the optical-electrical conversion coefficient rj can achieve 
higher ergodic capacity of MIMO FSO systems over i.n.i.d. 
rr turbulence channels. 


IV. Numerical Results 

In this section, the analytical results of MIMO FSO systems 
over different system conhguration scenarios and/or various 
channel fading conditions are presented and compared with 
Monte-Carlo simulations. The impact of some system and 
channel parameters on the ergodic capacity is also analyzed in 
detail. Moreover, the ergodic capacity for SISO FSO systems 
is also provided as a benchmark for comparison. We assume 
a wavelength of A = 850 nm and an optical-to-electrical 
conversion coefficient p = 1 without loss of generality. In 
Monte-Carlo simulations, 10^ i.i.d. and/or i.n.i.d. FT random 
samples are generated. In all the cases considered here, there 
is a good match between the analytical and the respective 
simulated results, which validates the accuracy of the proposed 
expressions. 

Figure [T] portrays the simulated, analytical (l22l l and high- 
SNR (l3Tt approximated ergodic capacity of MIMO FSO 
systems as a function of the average electrical SNR 7 for 
various values of M and N. We assume the strong turbulence 
fading channels of the MIMO FSO system are i.i.d., with 
hx values of = 3 x D = 0.01 m, and 

L = 4000 m. For comparison, we also plot the ergodic 
capacity of MIMO FSO systems over AWGN channel, whose 
ideal capacity is Cawgn = log 2 (1 + SNR), as the benchmark 
for comparison. It is clear to see that the analytical ergodic 
capacity is very accurate in the entire SNR regime, and 
the high-SNR approximation curves are quite tight in the 
moderate and high SNRs. As expected, the ergodic capacity is 
signihcantly improved as the number of transmit and receive 
apertures increases, especially in high-SNR regime. Finally, 
the gap between the MIMO AWGN and analytical curves 
decreases as M and N get larger. When M = 2 and A^ = 4, 
the analytical ergodic capacity is very close to the capacity of 
the MIMO AWGN channel. 

Recall that the MIMO FSO systems favor high data rate 
with broad spectrum, it is signihcant to achieve an accurate 



Eig. 1. Simulated, analytical, and high-SNR approximation ergodic capacity 
of MIMO ESO systems over i.i.d. distributed strong turbulence fading chan¬ 
nels versus the transmit SNR p (C^ = 3 x 10~^'^m“^/®, D = 0.01 m, 
L = 4000 m, A = 850 nm and r] = 1). Moreover, the capacity of the 
non-turbulent channel, e.g., MIMO AWGN channel, is plotted to provide a 
benchmark for comparison. 


ergodic capacity approximation to guide the practical system 
design. In order to emphasize the accuracy of our proposed 
results, the absolute difference between the analytical and 
simulated ergodic capacity is presented in Table |I] Moreover, 
the approximation errors of the ergodic capacity deduced in 
llSn are also presented for comparison. From Table |I] it is clear 
to see that our approximation outperforms the results proposed 
in ED in terms of approximation error. Considering the typical 
case ofM = W = 2, 7 = 1 dB, our analytical result is 
near hundredfold accurate compared with II 2 TII . Moreover, the 
approximation error tends to ascend as the average electrical 
SNR increases when using the method proposed in ED. while 
the approximation error of our a-p approximation is stable and 
small in the entire SNR regime. 

Figure |2] illustrates the analytical (l30l) . high-SNR dTSl) . 
simulated, and AWGN capacity of MIMO FSO systems over 
i.n.i.d. strong turbulence fading channels. The parameters fc; 
and mi is given by Q and (| 6 ]l, respectively, and the expected 
fading coefficient O; is randomly chosen subject to a real 
normal distribution with mean one and standard variance 0 . 1 . 
As it is clearly illustrated, the analytical results are accurate in 
the entire SNR regime, indicating the validity of our results. 
Moreover, the high-SNR expressions become sufficiently tight 
even at moderate SNRs (e.g., 5 dB) when deploying only two 
apertures at both the transmitter and receiver sides. Note that 
the gap between the AWGN and analytical curves becomes 
stable as M and N are moderately large. 

After validating the accuracy of our derived expressions, 
we then investigate the impact of system parameters on the 
ergodic capacity. Due to the space limit, we only consider 
the i.i.d. FF fading channels in Figs. [DID Note that similar 
hndings can be observed for the i.n.i.d. case. In Fig. [D the 
ergodic capacity of MIMO FSO systems under i.i.d. FF fading 
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TABLE I 

Performance Error Comparison BETWEEN Our Results and the Results Presented in in bit/s/Hz (C^ = 3 x D = 0.01 

M, L = 4000 M, A = 850 NM and t] = 1) 


7 (dB) 

Our Results 

Results in ETl 

M= 1,N = 2 

II 

to 

II 

to 

M = 2,N = 4 

M = 1,N = 2 

II 

to 

II 

to 

M = 2,N = A 

-5 

I.OOE-032 

5.00E-04 

5.00E-04 

4.70E-03 

2.22E-02 

2.61E-02 

-2 

8.00E-04 

8.00E-04 

1.30E-03 

1.21E-02 

3.12E-02 

2.82E-02 

1 

l.OOE-04 

3.00E-04 

7.00E-04 

2.13E-02 

3.75E-02 

3.04E-02 

4 

l.lOE-03 

2.10E-03 

9.00E-04 

3.11E-02 

4.40E-02 

3.11E-02 

7 

6.00E-04 

2.30E-03 

8.00E-04 

3.83E-02 

4.71E-02 

3.17E-02 

10 

8.00E-04 

1.60E-03 

1.20E-03 

4.47E-02 

4.81E-02 

3.15E-02 


* The scientific E-notation is used here, for example, l.OOE-03 is equal to 1 X 10 



Fig. 2. Simulated, analytical, and high-SNR approximation ergodic capacity 
of MIMO FSO systems over i.n.i.d. distributed strong turbulence fading 
channels versus the transmit SNR p (C^ = 3 X 10“^^m“^/®, D = 0.01 
m, L = 4000 m, A = 850 nm and r) = 1). Moreover, the capacity of the 
non-turbulent channel, e.g., MIMO AWGN channel, is plotted to provide a 
benchmark for compailson. 


conditions is illustrated for various values of the strength of 
the atmospheric turbulence The arrow in Fig. [3 means 
that the corresponding curves denote the value of from 
3 X to 1 X The similar method of 

arrows are used in Figs. |4] and | 6 ] In all cases, the link distance 
is hxed as L = 4000 m with constant D = 0.01 m and 
M = N = 2. The influence of the turbulence strength on the 
ergodic capacity becomes more obvious from weak to strong 
turbulence channels. For example, at 7 = 10 dB, the difference 
of the ergodic capacity between = 9 x 10“^^m“^/^ and 
— 1 X 10 “is much larger than that between 
Cl = i X 10“i%“2/3 and = 9 x lO-^^m-^/s. Fig. 
|3] indicates that the MIMO FSO systems almost achieve 
the AWGN capacity under clear weather or weak turbulence 
conditions (small values of Cl), while it losses much capacity 
in adverse weather or strong turbulence conditions (large 
values of Cl). 



Fig. 3. Simulated and analytical ergodic capacity of MIMO FSO systems 
over weak to strong i.i.d. distributed turbulence fading channels versus the 
transmit SNR p (M = N = 2, D = 0.01 m, L = 4000 m, A = 850 nm 
and T} = 1). Moreover, the capacity of the non-turbulent channel, e.g., MIMO 
AWGN channel, is plotted to provide a benchmark for comparison. 

The impact of the receiver aperture diameter D on the 
ergodic capacity of MIMO FSO systems over i.i.d. distributed 
strong turbulence fading channels is studied in Fig. |4] In 
addition, the AWGN capacity is plotted for comparison. The 
case of M = N = 2 apertures and L = 4000 m are 
considered. The graph likewise indicates that the analytical 
ergodic capacity coincides with simulated results in all cases 
under consideration, thereby validating the correctness of the 
our analytical results. Moreover, the receiver aperture diameter 
D has an obvious effect on the ergodic capacity, while its 
impact becomes more obvious as D becomes smaller. This 
observation can be explained by the fact that larger D tends to 
drive the fading channels approaching the deterministic case. 

In Fig. |5l the effects of optical link distance L on the 
ergodic capacity and scintillation index of MIMO FSO systems 
over i.i.d. distributed strong turbulence fading channels are 
investigated, assuming 7 = —5 dB, M = N = 2, D = 0.01 
m and Cl = fix 10“^‘*m“^/^. It is interesting to see that as L 
































Fig. 4. Simulated and analytical ergodic capacity of MIMO FSO systems and 
different receiver aperture diameter D over i.i.d. distributed strong turbulence 
fading channels versus the transmit SNR p {M = N = 2, = 3 x 

L = 4000 m, A = 850 nm and r] = 1). Moreover, the 
capacity of the non-turbulent channel, e.g., MIMO AWGN channel, is plotted 
to provide a benchmark for comparison. 



Distance of link L (m) 



Fig. 5. Simulated, analytical ergodic capacity and scintillation index of 
MIMO FSO systems over i.i.d. distributed strong turbulence fading channels 
versus the distance of link L (7 = —5 dB, M = N = 2, D = 0.01 m, 
7 = 3 X A = 850 nm and )? = 1). 


increases, the scintillation index in (|7]i approaches a maximum 
value greater than 1, and then changes the slope at a distance of 
about 4000 m. We observe that as scintillation index increases, 
the turbulence effect is getting stronger and thus the ergodic 
capacity decreases. Indeed, it is not a proportional relationship 
between the ergodic capacity and link distance as expected. 
This conclusion is consistent with the results presented in HU, 
EH. Therefore, the FSO system can achieve a certain capacity 
gain by adjusting the optical link distance. 

Assuming that the received optical signal undergoes i.n.i.d. 
turbulence channels, and the fading coefficient ili varies in 



Fig. 6. Simulated and analytical ergodic capacity of MIMO FSO systems and 
different fading coefficient ratio /3 over i.n.i.d. distributed strong turbulence 
fading channels versus the transmit SNR p {M = N = 2, = 3 x 

D = 0.01 m, L = 4000 m, A = 850 nm and r] = 1). 
Moreover, the capacity of the non-turbulent channel, e.g., AWGN channel, is 
plotted to provide a benchmark for comparison. 

each link. We consider a typical case of fl; = {2 < I < 

L), where /3 > 1 is the ratio of adjacent fading coefficients, 
to investigate the effect of fading coefficient volatility on 
the ergodic capacity. For the considered MIMO FSO system 
with AWGN receiver noise and EGC receiver, the relationship 
between the ergodic capacity and (3 is studied in Fig. |6] where 
the analytical results are generated from (1221) for i.i.d. FF 
channels (e.g., /3 = 1) and ( |30] | for i.n.i.d. FF channels (e.g., 
/3 > 1). For easy tractability but without loss of generality, 
all 0;s are normalized to yield an identical electrical SNR 
given by Q. From Fig. |6] it is clear that the ergodic capacity 
is descending as /3 increases (e.g., the ergodic capacity for 
/3 = 2 is much lower than that for /3 = 1), or equivalently, the 
ergodic capacity is shown to be larger when Oj from different 
links are more balanced, i.e., with lower volatility. Moreover, 
there is a significant reduction in ergodic capacity when one of 
the receive apertures collect most of the transmission power, 
such as /3 = 2 where 12^ > observation can 

be easily interpreted due to the EGC scheme used in our 
analysis. Finally, it is worthy to note from Fig. |6] that the 
analytical results are not very accurate with f] = 1.5 and 
/3 = 2. However, this approximation acts as a lower bound 
for all cases. 

V. Conclusions 

In this paper, the ergodic capacity performance of MIMO 
FSO systems with the EGC scheme turbulence channels was 
studied. An a-p approximation has been used for the sum of 
i.i.d. FF RVs, while a nested finite weighted approximation 
was employed for the sum of i.n.i.d. FF RVs. Analytical 
and accurate expressions of the ergodic capacity in forms of 
Meijer’s G- and Fox’s iF-functions were derived for both i.i.d. 
and i.n.i.d. cases. Accurate asymptotic ergodic capacity ex- 
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pressions at the high-SNR regime were also derived to provide 
useful insights regarding the parameters that affect the system 
performance. For example, a capacity gain can be achieved by 
using more transmit and receive apertures as well as increasing 
the diameter of receive apertures. Interestingly, for practical 
MIMO FSO systems, it is not a proportional relationship 
between the ergodic capacity and link distance as expected. 
Moreover, analytical results accompanied with Monte-Carlo 
simulations were presented to demonstrate effectiveness of the 
proposed expressions. Finally, the analysis method proposed 
in this paper can be used to derive the ergodic capacity of 
MIMO FSO system for other receiver scheme, such as MRC. 
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